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This  numerical  model  study  was  authorized  by  the  US  Army  Engineer  Dis¬ 
trict,  Portland  (NPP) ,  and  conducted  at  the  Coastal  Engineering  Research  Cen¬ 
ter  (CERC)  of  the  US  Army  Engineer  Waterways  Experiment  Station  (WES) .  The 
study  was  conducted  and  this  report  was  prepared  during  the  period  February- 
August  1985  by  Mrs.  Mary  A.  Cialone,  Research  Division,  CERC,  under  the  super¬ 
vision  of  Mr.  H.  Lee  Butler,  former  Chief,  Coastal  Processes  Branch,  present 
Chief,  Research  Division.  General  supervision  at  the  time  of  the  study  was 
provided  by  Dr.  R.  W.  Whalin,  former  Chief,  CERC,  Mr.  C.  C.  Calhoun,  Jr., 
former  Acting  Chief,  present  Assistant  Chief,  CERC,  and  Dr.  J.  R.  Houston, 
former  Chief,  Research  Division,  present  Chief,  CERC. 

The  main  purpose  of  this  study  was  to  model  the  coastal  current  regime 
in  the  vicinity  of  the  Yaquina  Bay  entrance  and,  in  so  doing,  determine  the 
velocities  of  the  tidal  and  wave-induced  currents  at  the  tip  of  the  north 
jetty.  As  part  of  this  study,  the  current  field  for  five  rehabilitative 
alternatives  was  modeled  to  provide  design  guidance  to  NPP  for  the  rehabili¬ 
tation  of  the  north  jetty  at  Yaquina  Bay. 

This  report  describes  the  application  of  three  numerical  models  to  the 
Yaquina  Bay  area.  The  numerical  computations  associated  with  this  work  were 
performed  on  the  CYBER  865  and  CYBER  205  Control  Data  Corporation  Cyber  com¬ 
puters  at  Rockville,  Maryland. 

Director  of  WES  during  this  investigation  and  during  the  preparation  and 
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publication  of  this  report  was  COL  Allen  r.  Grum,  USA.  COL  Dwayne  G.  Lee, 
CE,  is  the  present  Commander  and  Director.  The  Technical  Director  is 
Dr.  Robert  W.  Whalin. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Mul tipi 


degree  (angle) 


miles  (nautical  miles) 
miles  (US  statute) 


0.01745329 

0.3048 

1.852 

1.60934 


To  Obtain 


radians 
metres 
kilometres 
ki lometres 
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YAQUINA  BAY,  OREGON,  TIDAL  AND  WAVE-INDUCED 
CURRENTS  NEAR  THE  JETTIED  INLET 


Numerical  Model  Investigation 


1 .  Yaquina  Bay  is  located  on  the  northern  part  of  the  Oregon  coast 
about  110  miles*  south  of  the  Columbia  River  (Figure  1).  The  harbor  entrance 
is  protected  by  two  rubble-mound  jetties  1000  ft  apart.  The  navigation  channel 
is  400  ft  wide  at  its  outermost  reach,  300  ft  wide  at  sta  0+00  (Figure  2)  and 
900  to  1,200  ft  wide  in  the  turning  basin.  The  entrance  channel  is  dredged  to 
40  ft  at  the  outermost  portion,  thence  tapering  to  a  depth  of  30  ft  at  sta  0+00. 

2.  The  entrance  is  flanked  by  Yaquina  Reef  to  the  north,  approximately 
4,000  ft  offshore.  The  narrow  rock  reef  extends  to  the  north,  approximately 
parallel  to  the  coast  for  6,000  ft.  To  the  south  is  a  less  prominent  feature 
called  South  Reef.  The  length  of  its  shallowest  part  is  about  1,650  ft,  ex¬ 
tending  southwest. 


The  Problem 


3.  The  rubble-mound  jetties  were  originally  constructed  in  1918  to 
improve  navigation  into  Yaquina  Bay,  with  the  entrance  at  a  natural  opening  in 
the  reef  system.  In  1958,  construction  was  initiated  to  extend  the  jetties  to 
their  present  length,  having  the  north  jetty  nearly  reaching  the  southern  tip 
of  Yaquina  Reef.  Construction  was  completed  in  1966,  but  by  1970  the  outer 
330  ft  of  the  north  jetty  had  collapsed.  Rehabilitation  of  the  north  jetty 
was  completed  in  1978.  In  1979  60  ft  had  again  collapsed,  in  1980  250  ft  had 
collapsed,  and  in  1984  360  ft  had  collapsed  below  mean  lower  low  water  (mllw) . 
Some  of  the  rubble  material  has  slumped  inward  (or  to  the  south)  of  the  north 
jetty.  The  current  between  Yaquina  Reef  and  the  north  jetty  is  believed  to  be 

*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 
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Figure  2.  Location  of  jettied  entrance  to  Yaquina  Bay 


* 


a  major  contributing  factor  in  the  undermining  and  subsequent  collapse  of  the 
jetty  tip.  The  problem  faced  by  the  US  Army  Engineer  District,  Portland 
(NPP) ,  is  designing  rehabilitative  alternatives  for  the  north  jetty.  The 
Coastal  Engineering  Research  Center  (CERC)  was  requested  by  NPP  to  provide 
guidance  in  the  jetty  design.  Specifically,  the  problem  faced  by  CERC  was  to 
define  the  local  current  regime  using  hydrodynamic  numerical  models. 
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Model  Approach 


4.  The  objective  of  this  study  was  to  describe  the  currents  in  and 
around  the  entrance  to  Yaquina  Bay  by  using  two  hydrodynamic  numerical  models: 

a.  The  numerical  tidal  circulation  model  used  is  the  US  Army  Engi¬ 
neer  Waterways  Experiment  Station's  (WES)  Implicit  Flooding 
Model  (WIFM)  (Butler  1980).  WIFM  solves  the  vertically  inte¬ 
grated  Navier-Stokes  equations  by  using  an  alterr.n  ing  direction 
implicit  (ADI)  finite  difference  scheme.  The  code  uses  a  vari¬ 
ably  spaced  finite  difference  grid  to  allow  detailed  modeling  of 
areas  that  are  important  to  the  physical  problem.  Small  fea¬ 
tures,  such  as  jetties,  are  approximated  as  thin  barriers  be¬ 
tween  computational  cells.  The  code  allows  for  flooding  and 
drying  of  low-lying  areas  as  the  surface  elevation  changes. 

b.  The  wave-induced  current  model  used  is  CURRENT.  This  model  was 
developed  at  WES  by  Dr.  S.  Rao  Vemulakonda  (1984).  CURRENT 
solves  the  vertically  integrated  Navier-Stokes  equations  by 
using  an  ADI  finite  difference  scheme  and  the  Longuet-Higgins 
approach  to  radiation  stress.  Mixing  and  advection  terms  are 
included.  The  code  employs  a  variably  spaced  finite  difference 
grid  to  allow  a  detailed  modeling  of  areas  that  are  important  to 
the  physical  problem.  CURRENT  can  approximate  small  features, 
such  as  jetties,  as  thin  barriers  between  computational  cells. 

5.  CURRENT  requires  the  wave  number,  wave  angle,  and  wave  height  at 
each  grid  cell  as  input.  Therefore,  the  Regional  Coastal  Processes  Wave 
Propagation  Model  (RCPWAVE)  is  run  for  each  wave  condition  required  for  the 
study.  RCPWAVE  predicts  the  transformation  of  monochromatic  waves  over  com¬ 
plex  bathymetry,  including  refractive  and  diffractive  effects  (Ebersole, 
Cialone,  and  Prater  1985). 

6.  CURRENT  is  run  until  a  steady-state  condition  is  reached.  Calibra¬ 
tion  and  verification  of  the  tidal  model  are  accomplished  by  adjusting  model 
parameters  until  computed  velocities  correspond  to  observed  velocities  at 
select  locations.  After  verification,  several  rehabilitative  alternatives  are 
modeled  to  determine  their  influence  on  the  current  regime. 
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PART  II:  THE  NUMERICAL  TIDAL  CIRCULATION  MODEL 


The  Model 

7.  The  numerical  tidal  circulation  model  employed  is  the  WES  Implicit 
Flooding  Model  (WIFM) .  The  Navier-Stokes  equations  in  a  Cartesian  coordinate 
system  are  used  to  derive  the  hydrodynamic  equations  employed  by  WIFM 
(Figure  3).  The  two-dimensional  (2-D)  equations  of  continuity  ar.d  momentum 


contained  ,in  WIFM  are  obtained  by  assuming  vertical  accelerations  are  small 
and  the  fluid  is  homogeneous: 

Continuity 


9V  ,  3v  3v 

—  +  u  —  +  v  — 

at  ax  av 


fu  +  g  (n  -  na)  +  (u 

c  d 


2  4.  2 

+  V 


,1/2 


( 3) 


These  equations  are  solved  for  n  ,  u  ,  and  v  ,  which  denote  the 
water  surface  elevation  and  the  vertically  averaged  velocities  in  the 

x-  and  y-directions ,  respectively.  The  independent  variables  include: 
t  =  time 

d  =  p  -  h  =  total  water  depth 
h  =  still-water  elevation 

R  -  rate  at  which  the  water  volume  in  the  system  changes  through 
rainfall  and  evaporation 

f  =  Coriolis  parameter 

g  =  acceleration  due  to  gravity 

n  =  hydrostatic  water  elevations  due  to  atmospheric  pressure 
differences 

c  =  Chezv  coefficient  for  bottom  friction 
e  =  eddy  viscosity  coefficient 

F  ,  F  =  external  forcing  functions  in  the  x-  and  y-directions, 

X  y  respectively  (i.e.,  wind  stress) 

8.  WIFM  employs  a  variably  spaced  finite  difference  grid  in  its  compu¬ 
tational  scheme.  The  grid  axes  are  aligned  such  that  the  y-axis  is  approxi¬ 
mately  parallel  to  the  shoreline,  increasing  in  the  southerly  direction,  and 
the  x-axis  is  approximately  perpendicular  to  the  shoreline,  increasing  in  the 
offshore  direction  at  Yaquina  Bay  (Figure  4).  The  grid  cells  range  in  size 
from  145  to  1,460  ft  in  the  x-direction  and  from  108  to  1,480  ft  in  the 
y-direction.  The  depth  at  each  of  the  55  by  49  grid  cells  is  digitized  from 
the  1982  Yaquina  Bay  and  River  Bathymetry  Survey  to  a  datum  of  mllw.  The  max¬ 
imum  depth  of  the  Yaquina  Bay  WIFM  grid  is  155  ft  below  mllw.  All  data  are 
then  converted  to  the  National  Geodetic  Vertical  Datum  (NGVD)  by  subtracting 
4.03  ft. 

9.  Boundary  conditions  specified  as  input  to  WIFM  may  include  open 
boundaries,  land-water  boundaries,  and  subgrid  barriers.  Open  boundary 
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Figure  4.  WIFM  Grid 

conditions  exist  when  the  edge  of  the  computation?!  grid  is  defined  as  water; 
in  the  case  of  a  seaward  boundary  or  when  a  channel  lies  on  the  edge  of  a  com¬ 
putational  grid,  WIFM  requires  that  the  water  elevation  or  velocity,  in  tabu¬ 
lar  or  constituent  form,  be  specified  as  forcing  functions  on  such  boundaries. 

10.  Land-water  boundaries  (i.e.,  the  coastline)  can  be  either  fixed  or 
variable  in  the  WIFM  code.  For  a  fixed  boundary,  a  no-flow  condition  (u  ■  0 
or  v  «  0)  is  specified  at  the  land-water  interface.  A  variable  land-water 
boundary  condition  allows  for  flooding  and  drying  of  low-lying  areas.  This  is 
accomplished  by  comparing  the  water  elevation  in  a  water  cell  with  the  land 
elevation  in  an  adjacent  "dry"  cell.  When  the  water  level  exceeds  the  land 
elevation,  water  is  moved  into  the  dry  cell  by  using  a  broad-crested  weir 


formula  (Reid  and  Bodlne  1968) .  After  a  specified  amount  of  water  has  inun¬ 
dated  the  dry  cell,  it  is  thereafter  treated  as  an  open  boundary  face  and  com¬ 
putations  for  n  »  u  ,  and  v  are  made.  The  process  of  drying  a  water  cell 
is  modeled  in  a  manner  opposite  to  the  process  just  described. 

11.  Subgrid  barriers  are  thin  walls  defined  along  cell  faces.  There 
are  three  types  of  thin-wall  barriers  modeled  by  WIFU:  exposed,  submerged, 
and  overtopping.  Exposed  barriers  require  no  flow  across  the  specified  cell 
faces.  Submerged  barriers  are  modeled  by  controlling  flow  across  a  cell  face 
by  using  a  time-dependent  friction  coefficient.  The  third  type  of  subgrid 
barrier  is  the  overtopping  barrier.  These  barriers  are  dynamic  (i.e.,  they 
can  become  either  completely  exposed,  completely  submerged,  or  they  can  act  as 
broad-crested  weirs  as  model  time  progresses).  The  type  of  subgrid  barrier 
used  in  the  model  is  a  function  of  its  height  and  the  water  elevation  in  the 
two  cells  adjacent  to  the  barrier.  The  WIFM  grid  boundaries  and  barriers  are 
displayed  in  Figure  5. 


Figure  5.  WIFM  grid  boundaries  and  barriers. 
The  shaded  regions  are  land  cells 


Prototype  Data 

12.  To  calibrate  and  verify  WIFM  for  tidal  events,  prototype  tidal  ele¬ 
vation  and  water  current  data  are  required.  Tidal  elevations  are  needed  at 
the  open  ocean  and  inside  the  bay  entrance  as  boundary  conditions  for  WIFM. 

13.  Oregon  State  University  (OSU)  supplied  the  established  tidal  datums 
at  South  Beach,  Oreg. ,  inside  the  bay  entrance,  and  hourly  water  elevations 
during  the  month  of  September  1984  at  the  same  location. 

14.  CERC  placed  a  Temperature-Depth  Recorder  (TDR)  near  the  70-ft  sea 
buoy  located  at  44°36'  N  latitude,  124°06'  W  longitude  from  24-26  September 
1984  for  48  hr.  The  tidal  elevations  were  recorded  relative  to  an  arbitrary 
datum  every  7.5  min.  The  OSU-established  datums  and  short-term  data  were  used 
to  correlate  the  CERC  data  to  the  NGVD  for  use  as  an  open-ocean  boundary  con¬ 
dition  in  the  tidal  model. 

15.  The  US  Department  of  Commerce  1984  provided  tabulated  tidal  con¬ 
stituents  at  South  Beach  (44°38'  N  latitude,  124°03'  W  longitude)  that  could  be 
used  by  WIFM  as  a  channel-boundary  condition.  The  eight  constituents  used  for 
the  tide  at  South  Beach  (M2,  Kl,  01,  S2 ,  N2,  PI,  K2,  and  Ql)  form  90  percent  of 
the  total  tide.  Table  1  provides  the  amplitudes  and  epochs  for  the  eight  con¬ 
stituents  aforementioned. 

Table  1 

Tidal  Constituents  at  South  Beach,  Oreg. 


Constituent 

Height,  ft 

Kelvin,  deg 

M2 

2.884 

342.97 

Kl 

1.431 

113.83 

01 

0.878 

96.74 

S2 

0.  782 

10.52 

N2 

0.595 

319.96 

PI 

0.439 

111.83 

K2 

0.230 

12.90 

Ql 

0.161 

88.95 

12 


*  i. 


16.  Current  data  were  taken  in  the  inlet  along  two  ranges  (six  gage 
locations)  using  an  Endeco  Model  110  current  meter  (Figure  2).  Hourly  read¬ 
ings  were  recorded  by  CERC  during  an  11-hr  period  on  25  September  1984.  A 
second  day's  field  data  collection  effort  on  26  September  1984,  at  the  six 
gage  locations,  proved  to  be  unsuccessful.  Therefore,  a  limited  amount  of 
data  were  available  for  calibration  of  the  tidal  model . 

Calibration 

17.  Calibration  of  WIFM  begins  with  the  selection  of  boundary  condi¬ 
tions  for  all  open-water  cells  at  the  edges  of  the  computational  grid.  In 
this  study,  tidal  elevations  were  used  as  forcing  functions  on  the  seaward 
boundaries  and  the  inlet  boundary.  WIFM  is  calibrated  by  adjusting  the  bound¬ 
ary  conditions,  the  Chezy  frictional  coefficients  in  each  cell,  and/or  the 
bathymetry  until  the  model  computes  the  correct  tidal  signals  or  velocities 
for  locations  inside  the  grid.  The  model  was  required  to  reproduce  the  veloc¬ 
ities  at  six  gages  in  the  inlet  during  an  11-hr  period  on  25  September  1984 
for  this  study.  Since  the  field  data  collection  effort  was  not  extensive, 
this  single  day  of  data  was  used  for  calibration  and  verification  of  WIFM. 

18.  The  calibration  procedure  started  by  using  the  TDR  prototype  tidal 
elevations  on  the  open-ocean  boundaries  and  modified  TDR  prototype  tidal  ele¬ 
vations  at  the  inlet  boundary.  The  inlet  boundary  condition  was  modified  be¬ 
cause  its  tidal  signal  is  expected  to  lag  that  of  the  seaward  boundary.  The 
frictional  values  in  the  inlet  were  adjusted  to  account  for  shoals  and  the 
jetty  material  which  slumped  into  the  inlet.  With  these  adjustments,  the 
numerical  model  tidal  velocity  predictions  adequately  matched  the  prototype 
tidal  velocities  (Plates  1-3). 

19.  Discrepancies  occurring  at  sta  A  and  sta  2  may  be  caused  by  the 
lack  of  data  concerning  the  amount  of  slumping  near  the  north  jetty.  Another 
factor  which  may  have  contributed  to  the  discrepancy  is  the  cell  size  in  this 
region.  Although  fine  cells  (150  ft)  were  used  in  the  north  jetty  region,  the 
grid  may  not  have  been  fine  enough  to  resolve  the  velocities  accurately. 
Running  WIFM  with  smaller  grid  cells  was  economically  impractical  due  to  the 
small  time-step  the  grid  would  require. 


PART  III:  THE  RECIONAL  COASTAL  PROCESSES 
WAVE  PROPAGATION  MODEL 


The  Model 

20.  The  RCPWAVE  model  predicts  the  transformation  of  monochromatic 
waves  over  complerx  bathymetry,  including  refractive  and  diffractive  effects. 
The  wave  propagation  model  computes  the  wave  height,  direction,  and  number 
throughout  a  finite  difference  grid.  These  wave  parameters  are  then  used  as 
forcing  functions  to  drive  CURRENT. 

21.  Berkhoff  (1972)  derived  an  elliptical  equation  to  approximate  the 
wave  transformation  of  linear  waves  over  an  arbitrary  bathymetry  with  mild 
bottom  slope.  The  governing  equation  is  given  by 


where 

c(x,y)  =  wave  celerity 

c  (x,y)  =  group  velocity 

®  Is 

<Kx,y)  =  the  complex  velocity  potential  =  ae 

a(x,y)  =  a  wave  amplitude  function 

e  =  2.71828 

i  =  vCT 

s(x,y)  =  wave  phase  function 

o  =  angular  wave  frequency 

The  wave  number  is  found  from  the  dispersion  relation 

=  gk  tanh  (kh) 

where 

g  =  acceleration  due  to  gravity 
k(x,y)  =  wave  number 
h(x,y)  «  still-water  depth 

By  substituting  the  expression  for  the  velocity  potential  into  Equation  4  and 
solving  for  the  real  and  imaginary  parts,  two  equations  can  be  derived: 
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These  equations,  along  with  the  dispersion  relation  and  the  irrotationality 
condition,  describe  the  combined  refraction  and  diffraction  process. 

22.  The  RCPWAVE  model  has  the  capability  of  computing  solutions  on  rec¬ 
tilinear  grid  systems  having  uniform  or  variably  sized  grid  cells.  The  coor¬ 
dinate  system  and  grid  cell  convention  used  in  the  model  are  given  in  Fig¬ 
ure  6.  Lateral  boundary  conditions  are  such  that  the  value  at  cells  J  =  N 
and  J  =  1  are  set  equal  to  the  value  at  cells  J  =  N  -  1  and  J  =  2  , 
respectively.  Boundary  conditions  along  the  seaward  length  of  the  grid  are 
computed  from  deepwater  wave  conditions  and  are  used  to  initiate  the  shoreward 
marching  scheme  of  the  model,  assuming  Snell's  law  is  valid  from  deepwater  to 


the  outer  boundary  of  the  grid. 
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Figure  6.  Definition  of  coordinate  system  and  grid  cell 
conventions  used  in  the  model 
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Model  Approach 


23.  The  RCPWAVE  model  is  run  on  two  separate  grids.  The  larger  (93  by 
83),  variably  spaced,  rectangular  grid  used  for  this  model  is  shown  in  Fig¬ 
ure  7.  This  grid  was  developed  for  the  wave  propagation  analysis  performed  by 
CERC  as  part  of  the  Yaquina  Bay  North  Jetty  Rehabilitation  Study.  Grid  dimen¬ 
sions  are  20.1  nautical  miles  (n.m.)  in  the  longshore  (y)  direction  and 
18.6  n.m.  in  the  offshore  (x)  direction,  with  the  y-axis  6.5  deg  east  of  true 
north  (Figure  8).  Water  depths  are  digitized  from  the  NOS  nautical  charts, 
except  near  the  jetties  where  NPP  provided  survey  data.  Depths  relative  to 
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Figure  7.  RCPWAVE  93  by  83  grid 
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Figure  8.  RCPWAVE  93  by  83  grid  boundaries 

mllw  range  fron  approximately  -12.0  to  -40.0  ft  along  the  y-axis  (I  *  1)  and 
-300.0  to  -500.0  ft  along  the  seaward  boundary  (I  •  93).  Deepwater  wave  in¬ 
formation  is  required  as  input  to  the  model,  specifically  height,  period,  and 
direction.  Table  2  displays  the  wave  conditions  selected.  The  model  is  run 
on  the  larger  grid  to  propagate  waves  from  far  offshore  into  the  jetty  area. 
For  each  of  the  24  runs,  the  wave  height  and  angle  at  each  grid  cell  are 
stored  permanently. 
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Deepwater  Wave 
Height,  ft 

14 
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Table  2 

RCPWAVE  Model  Wave  Conditions 


Period,  sec 

8.0 

11.0 

14.0 

16.0 

8.0 

11.0 

14.0 

16.0 


Direction,  deg 


West-northwest  (-16.0)* 


West  (6.5) 
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14.0 

16.0 

8.0 

11.0 

14.0 

16.0 

8.0 

11.0 

14.0 

16.0 

8.0 

11.0 

14.0 

16.0 


West-southwest  (29.0) 


West-northwest  (-16.0) 


West  (6.5) 


West-southwest  (0.0) 


24.  The  stored  values  of  wave  height  and  angle  are  then  Interpolated 
onto  a  finer,  smaller  grid  (Figure  9).  The  interpolation  program  uses  a 
fairly  simple  scheme  to  determine  values  of  wave  height  and  angle  at  each  grid 
cell  in  the  smaller,  "new"  grid.  For  every  new  grid  cell,  a  search  is  done  in 
each  90-deg  quadrant  for  the  nearest  value  of  wave  height  from  the  larger 
"old"  grid.  A  weighted  average  of  these  four  values  is  used  to  compute  a  new 
wave  height.  The  weighting  functions  are  based  on  relative  distances  from  the 
new  grid  centers  to  the  cell  centers  defining  the  positions  of  the  four 
closest,  old  wave  height  values.  The  procedure  is  repeated  to  determine 
values  of  the  wave  angle  throughout  the  new  grid. 

25.  The  RCPWAVE  model  is  run  on  the  smaller,  finer  grid.  The  initial 
guess  of  wave  heights  and  wave  angles  over  the  entire  grid  is  the  interpolated 
results  from  the  larger  grid  rather  than  the  results  from  using  Snell's  law. 
The  smaller  (61  by  80),  variably  spaced,  rectangular  grid  used  by  the  model  is 
shown  in  Figure  9.  Grid  dimensions  are  4.58  n.m.  in  the  longshore  (y)  direc¬ 
tion  and  2.05  n.m.  in  the  offshore  (x)  direction,  with  the  y-axis  9.5  deg  east 
of  true  north.  Water  depths  are  digitized  from  the  1982  Yaquina  Bay  and  River 
Bathymetry  Survey  except  in  the  entrance  channel  where  NPP  provided  a  more 
detailed  bathymetric  chart.  Depths  relative  to  mllw  range  from  approximately 
-1.0  to  -96.0  ft.  Since  the  RCPWAVE  model  does  not  model  structures,  a  rou¬ 
tine  was  developed  to  use  with  the  RCPWAVE  model  to  mode]  the  diffractive 
effect  of  the  jetties.  The  wave  model  is  run  on  the  smaller  grid  to  provide 
the  wave  height,  angle,  and  number  at  each  grid  cell  as  input  to  CURRENT. 
CURRENT  also  uses  the  61  by  80  cell  grid;  therefore,  no  interpolation  is 
required . 


PART  IV:  THE  WAVE-INDUCED  CURRENT  MODEL 


The  Model 


26.  CURRENT  is  used  to  predict  the  currents  induced  by  breaking  waves 
at  Yaquina  Bay.  The  model  takes  the  radiation  stress  approach  of  Longuet- 
Higgins  and  solves  the  vertically  integrated  equations  of  continuity  and 
momentum  by  using  an  ADI  scheme.  Mixing  and  advection  terms  are  included. 

The  hydrodynamic  equations  used  in  the  model  are  derived  from  the  Navier- 
Stokes  equations,  assuming  the  fluid  is  homogeneous  and  incompressible  and  the 
vertical  accelerations  are  negligible.  The  2-D,  depth-averaged  forms  of  the 
equations  are  obtained  by  vertically  integrating  the  three-dimensional  form  of 
the  equations  with  the  appropriate  boundary  conditions  applied.  The  2-D  equa¬ 
tions  are  time-averaged  over  one  wave  period  to  give 
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where 


u  and  v  =  depth-averaged  horizontal  velocity  components  at 
time  t  in  the  x-  and  y-directions ,  respectively 

g  =  acceleration  due  to  gravity 

n  =  displacement  of  the  mean-free  surface  with  respect 
to  the  still-water  level 


T,  and 
bx 


0 

d 

h 


Tby 


=  mass  density  of  seawater 

=  n  -  h  •  total  water  depth  (Figure  10) 

=  bed  elevation  with  respect  to  the  still-water  level 
((+)  land,  (-)  water) 

=  bottom  friction  stresses  in  the  x-  and  y-directions, 
respectively 


a.  CROSS  SECTION  A-A 
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Figure  10.  Definition  sketch  for  an  irregular  beach 
(swl  ”  still  water  level,  mwl  ■  mean  water  level) 

S  ,  and  S  ■  radiation  stresses  which  arise  because  of  the 
^  excess  momentum  flux  due  to  waves 

t  =  lateral  shear  stress  due  to  turbulent  mixing 
xy 

Initial  and  Boundary  Conditions 


27.  To  run  CURRENT,  appropriate  initial  and  boundary  conditions  must  be 
specified.  At  Yaquina  Bay,  the  model  calculations  are  started  from  rest, 
setting  n  ,  u  ,  and  v  to  zero.  Model  boundaries  are  shown  in  Figure  11. 
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Figure  11.  CURRENT  and  RCPWAVE  grid  boundaries 

The  radiation  stress  gradients  are  gradually  applied  to  avoid  unwanted  tran¬ 
sients  in  the  solution.  This  is  accomplished  by  gradually  building  up  to  the 
full  radiation  stress  values  over  a  number  of  time-steps.  The  numerical  com¬ 
putations  continue  until  a  steady-state  condition  is  reached. 

28.  CURRENT  permits  the  following  boundary  conditions: 

a.  "No  flow"  (wall) .  This  boundary  condition  is  used  at  closed 
boundaries  (i.e.,  land-water  boundaries  and  at  impermeable 
structures).  Land-water  boundaries  are  closed  because  CURRENT 
does  not  permit  flooding  and  drying.  The  jetties  at  Yaquina 
Bay  are  treated  as  impermeable,  nonovertopping  structures,  and 
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the  no-flow  boundary  condition  is  applied.  The  velocity  per¬ 
pendicular  to  closed  boundaries  is  set  to  zero. 

Elevation  (tide) .  This  condition  may  be  used  at  the  offshore 
boundary  since  the  setup  is  an  unknown  quantity.  It  should  be 
noted,  however,  that  a  radiation  stress  boundary  condition  is 
preferred  at  the  offshore  boundary. 

Uniform  flux.  The  uniform  flux  boundary  condition  is  used  for 
lateral  boundaries.  The  flux  at  a  boundary  cell  is  set  equal 
to  its  adjacent  interior  cell.  Thus,  the  condition  assumes 
3(ud)/3x  =  0  or  3(vd)/3y  =  0  at  the  boundary. 

Radiation.  The  radiation  boundary  condition  requires  that  any 
transients  initially  developing  within  the  numerical  grid 
propagate  out  of  the  grid  as  gravity  waves.  The  radiation 
condition  takes  the  following  form: 


+  c  in  -  o 

3t  3x 


where 

n(t,x)  =  surface  disturbance 
c  =  phase  speed  of  n 

CURRENT  uses  this  boundary  condition  at  the  offshore  boundary  rather  than 
boundary  condition  of  elevation  and  no  flow,  due  to  the  reflective  nature  of 
the  elevation  and  no-flow  boundary  conditions.  The  elevation  and  no-flow 
boundary  conditions  tend  to  trap  transients  between  the  offshore  and  nearshore 
boundaries,  and  damping  of  transients  would  take  place  at  a  slow  rate.  The 
radiation  boundary  condition  allows  the  transients  to  propagate  out  of  the 
grid,  producing  an  appropriate  n  value  at  the  offshore  boundary.  Figure  12 
displays  the  boundary  conditions  used  for  the  Yaquina  Bay  CURRENT  runs. 

Model  Approach 


29.  CURRENT  is  run  at  Yaquina  Bay  to  predict  the  currents  caused  by 
breaking  waves.  The  model  requires  the  wave  heights,  angles,  and  numbers 
throughout  a  gridded  region  as  input  conditions.  This  information  is  provided 
by  the  RCPWAVE  model.  The  24  deepwater  wave  conditions  used  to  run  RCPWAVE 
and  then  CURRENT  are  shown  in  Table  2. 

30.  CURRENT  employs  a  variably  spaced  finite  difference  grid  in  the 
computational  scheme.  The  grid  axes  are  aligned  such  that  the  y-axis  is  ap¬ 
proximately  parallel  to  the  shoreline,  increasing  in  the  southerly  direction 
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Figure  12.  CURRENT  model  boundary  conditions 

and  the  x-axis  is  approximately  perpendicular  to  the  shoreline,  increasing  in 
the  offshore  direction  at  Yaquina  Bay  (Figure  9).  The  RCPWAVE  model  was  also 
run  on  this  grid.  The  grid  cells  range  in  size  from  60  to  545  ft  in  the 
x-direction  and  from  150  to  511  ft  in  the  y-direction.  The  depth  at  each  of 
the  61  by  80  grid  cells  is  digitized  from  the  1982  Yaquina  Bay  and  River 
Bathymetry  Survey  to  a  datum  of  mllw.  The  maximum  depth  of  the  Yaquina  Bav 
CURRENT  grid  is  96  ft  below  mllw.  All  data  are  then  converted  to  NGVD  by  sub 
tracting  4.03  ft.  At  this  point,  CURRENT  is  run  until  a  steady-state  con¬ 
dition  is  reached. 


PART  V:  GRID  CONNECTION 


31.  The  purpose  of  this  study  is  to  model  the  velocities  of  the  tidal 
and  wave-induced  currents  in  the  vicinity  of  the  Yaquina  Bay  entrance  and,  in 
so  doing,  determine  the  current  regime  at  the  tip  of  the  north  jetty.  This  is 
accomplished  by  superimposing  the  tidal  and  wave-induced  currents  predicted  by 
the  hydrodynamic,  numerical  models,  WIFM  and  CURRENT.  Since  WIFM  is  run  on  a 
grid  of  49  by  55  cells  and  CURRENT  is  run  on  a  grid  of  61  by  80  cells,  inter¬ 
polation  is  required  before  superposition  can  be  performed.  The  49  by  55  WIFM 
grid  velocities  are  interpolated  to  the  finer  CURRENT  grid  by  using  the  fol¬ 
lowing  method: 

a.  For  every  (I,  J)  cell  in  the  CURRENT  grid,  a  search  is  executed 
for  the  nearest  WIFM  u-velocity  component  in  each  90-deg 
quadrant  of  the  WIFM  grid. 

b.  A  weighted  average  of  the  four  closest  WIFM  u-values  is  used 
to  compute  a  WIFM  u-value  for  the  CURRENT  (I,  J)  grid  cell. 
The  weighting  functions  are  based  on  the  relative  distances 
from  the  CURRENT  (I,  J)  grid  cell  to  the  cell  centers  defining 
the  positions  of  the  four  closest  WIFM  values. 

c.  The  procedure  is  repeated  to  determine  the  WIFM  v-velocity 
component  throughout  the  CURRENT  grid. 

32.  The  WIFM  values  of  u  interpolated  onto  the  CURRENT  grid  are 
superimposed  on  the  CURRENT  u-values  ,  providing  the  total  u-velocity  com¬ 
ponent  at  each  grid  cell.  The  v-velocity  is  dealt  with  in  the  same  manner. 
Vectorially  summing  u  and  v  ,  the  total  current  regime  at  Yaquina  Bay  is 
determined . 


26 


PART  VI:  REHABILITATIVE  ALTERNATIVES 


33.  Five  rehabilitative  alternatives  were  modeled  to  provide  the  cur¬ 
rent  field  throughout  the  gridded  region.  The  following  structural  alterna¬ 
tives  were  included: 

a.  Extending  the  north  jetty  to  its  original  length. 

b.  Adding  a  spur  near  the  end  of  the  north  jetty,  approximately 
perpendicular  to  Yaquina  Reef. 

c.  Adding  a  spur  near  the  end  of  the  north  jetty,  approximately 
perpendicular  to  the  jetty. 

d.  Extending  the  north  jetty  to  its  original  length  and  adding  a 
spur  near  the  end  of  the  north  jetty,  approximately  perpendicu 
lar  to  Yaquina  Reef. 

e.  Adding  an  L-shaped  barrier  near  the  end  of  the  north  jetty. 

34.  The  tidal  model  was  run  for  the  five  rehabilitative  alternatives. 
Structural  alternatives  a,  b,  d,  and  e  were  effective  in  altering  the  water 
current  pattern.  Alternatives  b  and  e  were  selected  to  use  with  CURRENT  be¬ 
cause  those  structural  alternatives  produced  beneficial  effects  on  the  wave- 
induced  current  pattern  near  the  north  jetty  tip. 


PART  VII:  RESULTS 


35.  After  calibration  and  verification,  the  tidal  model  was  run  with 
five  rehabilitative  alternatives  to  test  their  effects  on  the  tidal  current 
patterns  near  the  inlet.  The  structural  alternatives  were  modeled  as  thin- 
wall  barriers  in  the  tidal  model.  Plates  4  and  5  show  the  tidal  current  pat¬ 
terns  at  peak  flood  and  ebb  for  the  existing  conditions  at  Yaquina  Bay.  The 
tidal  current  patterns  at  peak  flood  and  ebb  for  the  five  rehabilitative  al¬ 
ternatives  are  displayed  in  Plates  6-15.  Plates  16-25  show  (vectorially)  the 
difference  between  the  existing  tidal  current  patterns  and  the  tidal  current 
patterns  produced  by  the  five  rehabilitative  alternatives  at  peak  flood  and 
ebb.  Structural  alternatives  a,  b,  d,  and  e  (Part  VI)  alter  the  tidal  current 
pattern  near  the  tip  of  the  north  jetty. 

36.  CURRENT  was  used  to  predict  the  existing  wave-induced  currents  at 
the  entrance  to  Yaquina  Bay  for  24  wave  conditions  (Table  2) .  Plates  26 
through  31*  show  the  steady-state,  wave-induced  current  patterns  produced  by 
CURRENT  for  6  of  the  24  wave  conditions.  Wave  period  did  not  appear  to  have  a 
strong  effect  on  the  wave-induced  current  patterns,  whereas  the  wave  height 
and  angle  influenced  the  current  patterns  significantly.  (Plates  are  shown 
only  for  the  various  wave  heights  and  angles.)  Plates  32-37  and  38-43  show 
the  superposition  of  tidal  and  wave-induced  currents  at  peak  flood  and  peak 
ebb,  respectively.  The  consistent  circulation  pattern  to  the  south  of  the 
south  jetty  is  related  to  the  acute  angle  that  the  structure  makes  with  the 
shoreline.  Coast  Guard  personnel  verified  the  existence  of  the  circulation 
pattern.**  The  wave-induced  currents  are  significantly  stronger  than  the 
tidal  currents  at  the  tip  of  the  north  jetty. 

37.  Testing  the  five  rehabilitative  alternatives  for  24  wave  conditions 
was  economically  infeasible.  Alternatives  b  and  e  were  selected  from  the  five 
rehabilitative  alternatives  to  use  with  CURRENT.  The  selection  was  based  on 
the  knowledge  of  the  wave-induced  current  patterns  produced  by  the  initial 

24  CURRENT  runs  for  existing  jetty  conditions  and  the  performance  of  the 

*  In  the  plates,  H  ■  deepwater  wave  height,  T  -  period,  0  =  wave 
direction. 

**  Personal  Communication,  December  3,  1985,  David  Edwards,  Boatswains  Mate 
First  Class,  US  Coast  Guard. 


five  rehabilitative  alternatives  in  varying  the  tidal  current  patterns.  Three 
wave  conditions  were  selected  for  testing  with  the  rehabilitative 
alternatives: 


Height 

Period 

Direction 

ft 

sec 

deg 

14 

11 

West-northwest 

14 

14 

West 

14 

11 

West-southwest 

The  three  wave  conditions  were  selected  to  represent  waves  arriving  from  the 
three  directions  and  occurrences  of  strong  currents  near  the  north  jetty. 

Using  a  wave  height  of  14  ft  best  produced  the  specified  conditions. 

Plates  44-61  display  the  wave-induced  current  patterns  and  the  superposition 
of  tidal  and  wave-induced  currents  at  peak  flood  and  peak  ebb  for  alternatives 
b  and  e. 

38.  In  comparing  Plates  32-34  with  Plates  47-49  and  56-58,  it  is  evi¬ 
dent  that  structural  alternatives  b  and  e  diminish  the  intensity  of  the  cur¬ 
rents  north  of  the  north  jetty  for  flood  tidal  flow.  At  ebb  tide,  a  compari¬ 
son  of  Plates  38-40  with  Plates  50-52  and  59-61  shows  a  similar  trend  in 
reducing  current  intensity  when  structural  alternatives  b  and  e  are  in  place. 
Currents  at  the  north  jetty  are  reduced  because  the  structures  redirect  the 
flow  outside  the  north  jetty  tip.  Flows  north  of  the  Inlet  are  generally 
north  to  south  and  parallel  to  the  Yaquina  Reef,  regardless  of  the  deepwater 
wave  direction.  The  current  appears  to  be  channelized  by  the  reef  and  struc¬ 
tural  alternatives  b  and  e  redirect  the  channelized  current. 


PART  VIII:  CONCLUSIONS 


39.  WIFM  provided  the  tidal  current  field  near  the  entrance  to  Yaquina 
Bay  for  existing  conditions  and  for  the  following  five  rehabilitative 
alternatives : 

a.  Extend  the  north  jetty  to  its  original  length. 

I 

b.  Add  a  spur  to  the  north  jetty,  approximately  perpendicular  to 
Yaquina  Reef. 

c.  Add  a  spur  to  the  north  jetty,  approximately  perpendicular  to 
the  north  jetty. 

d.  Extend  the  north  jetty  to  its  original  length  and  add  a  spur  to 
the  north  jetty,  approximately  perpendicular  to  Yaquina  Reef. 

e.  Add  an  L-shaped  barrier  to  the  north  jetty. 

The  results  of  this  model  investigation  indicate  that  rehabilitative  alterna¬ 
tives  a,  b,  d,  and  e  are  effective  in  altering  the  tidal  current  pattern  near 
the  tip  of  the  north  jetty. 

40.  CURRENT  provided  the  wave-induced  current  field  near  the  entrance 
to  Yaquina  Bay  for  24  wave  conditions.  For  each  wave  condition,  the  wave- 
induced  current  field  was  superimposed  onto  the  existing  tidal  current  field 
at  peak  flood  and  peak  ebb,  providing  the  total  current  field  at  Yaquina  Bay. 
CURRENT  was  then  run  for  three  selected  wave  conditions  with  rehabilitative 
alternatives  b  and  e.  For  the  selected  alternatives  b  and  e,  the  wave-induced 
current  field  was  superimposed  onto  the  tidal  current  field.  The  resulting 
current  patterns  demonstrate  that  structural  alternatives  b  and  e  are  effec¬ 
tive  in  altering  the  current  field  at  the  i ip  of  the  north  jetty. 
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Plate  6.  Tidal  current  pattern:  alternative  a  Part  VI;  peak  flood 
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Plate  10.  Tidal  current  pattern:  alternative  c;  peak  flood 


A4444A4AAAA  ►  > 

AA4A4AAAAA  4  >  > 

A44A4AAA  A  A  A  A  A 

A4444AAAAA  A  A  A 

AAA4AAAA  A  ►  A  A  a 

4A4AAAAA  A  A  A  T  T 

AA44AAAAAAtrT 
A4444AAA  A  A  A  A  « 


4  A  A  A 


A  A  A  A  A 
A  A  A  A  A 


AAA  A 


A  A  44A444444A4AAAAAA  A  A  A  A 
A  A  A444444444AAAAAAA  A  AAA 
A  *  # 444444444AAAAA  A  A  A  A  A  A 
A  A  *444444444«4«AA  A  A  A  A  A  A 
A  A  4  444 4444 4 ******  *  A  A  A  A  A 
A  4  ttf  #4444444******  A  A  A  A 

A  t  Tf  t  f  ff^44***««*  A  A  A  AAA 

T  •  •  a  « 

»>»■» ■»«  *  J 


A  A  >  A 


/  xn 


*  f-4T>  v 

*  A/.'. 


\/s 


fgr 

,x  x*  -»rr  *. 
x*  srr  ► 


„  x 

7*  x*  xV^£ 

^  x* 

X*  x'  x^> 

rM 


X  X 

/  . 


■3  *  *  V  iVwY***“*M* 

*  I  ^  «  UVliVWMMU 

v  «  U1IWM<<<<< 


mm 

P#T : 

Tr/fV' 

y  '•  . 

'  ^  0  A  *  >  > 


*>>>»>-*  4  A 

>>>>444  4 


>>>444  4 


>>>  4  4  4  4 


>>>>444  4 


>>>4444  4 


^M<mm«U444  4  4  4  4 

C  <<<<>>AAAAAAAAA A  A  4  4  4  4 
4  y  4  4  4  4  4 

4.  4.4.4.U'VrvAA*yy*4  4  A  A  4  4  4 
>  4>l.|pVtfvWWVAAA4  4  4  4  4  4 
U  i>t>|>t>VVvvvVW4V<  4  4  4  4  4  4 

V  VWWvwwvv  W44  4  4  4  4  4  4 
A  VVVWVVWVWV4  4  4  4  4  4  4  4 

V  VVVVWVV4444444  4  4  4  4  4  4 

V  VV44444V4V444444  4  4  4  4  4 
4  444444 44444444 4  4  4  4  4  4  4 
A  444444444444444  4  4  4  4  4  4 
4  444444444444444 4 4  4  4  4  4 
4  4444444444444444  4  4  4  4  4 
4  4444444444444444  4  4  4  4  4 


4  4  4444444444444444  4  4  4  4  4 


4  4  4  4444444444444444  4  4  4  4  4 


4  A  444444444444444  4  4  4  4  4  4 


4  4  4444444444444444  4  4  4  4  4 


■>?>?> 
v  , 

*  .  *  k> 


*.  \  •  v 


o.-x 

;  «  _  »T- 


grvJ>*K<fg<7^V  '  »  vj 


CvWWgEM'Vg1 


VWWVyvvvy 
V  |r  y  V  v  V  V  V  V 


Vvk-vvv  v  y 
y  y  W  y  v  v  y  y 


W  4  y  V  y  v  y  y  y 
t'fcWVyy  y  y  y 


y  y  y  y  y  y  v  y  y 

v  y  y  y  v  v  v  y  y 


vvwyyy  v  y  y 

vyyvyvvyv 


v  y  y  y  v  v  y  v  v 
vyyyyyvv  y 


vyyyyvwyy 

yi.4-i.yyyy 


v  v  4  4  4.  4. 

y  V  *  9  »»»»* **14.44.44  4  4  4  4  4  4 


4  4  4  4  4 

»» »»»»»»»4  4»444  4  4  4  4  4  4 


44  4  4  4  4  4 


•  ♦***MMi«**4 


w  y  y  y  y  y  y  y 

y  y  y  y  y  y  y 


y  y  y  w  44 

v  y  v  y  y 


y  y  y  4 


: 


7//yy^  '  a* 
Y/Y/  /  '  rr^ 

77/  /  "  -  -  - 

/y  /  /  r^J  4  A  *  *  * 

//  /  '  *  4  *  '  A  *  * 

/  t  *  *  «  «  *  *  *  * 


V/ 

^  /  s 
/  /  / 

sss 
ss  s 

/S' 
/  *  n- 


*1  <1  4  <4 


**WW***“4**4*  *  4  4  * 

■*  *  v  •« « « « . 

*  r—  l  [MU'**™*****  *  ♦  ♦  *  « 

I  *■  syj////M*****~~’*'''**  * 

S 

$yx%& . . 

y  ^✓<^r*44**m*****  * *  * 

y  ^  ^  |y  ^>^4  444444444  4  4  4  4  4  4 

.^Jy  Km-,AAAAaAAA444  4  A  4  4  4  4 
^  I  4  *•  rrr44AAA44444  4  4  4  4  4  4  4 

4  I  A4A4AAAAA444  4  4  4  4  4 

J  t  y  >  >>»>T't444444A44  4  4  4  4  4  4 
4  t  T  7>»,HmAAAAAA*A  4  4  4  4 
-f  ft«rrT*»-1i^^AAAA4  4  4  4  4  4  4 
*»  -»  I  5-tmmiAAAAAA»*  A  4  4  4 

A  44  11^111114AA4AAA  A  A  A  4  4  4 
A  4  A  AAAAAAAAAAAAAAAA4  A  4  4  4 
4  4  4  44444444AAAAAA  A  A  A  A  4  4  4 

A  4  4  444444444AAAAAA  A  A  A  4  A  A 

4  4  4  AAAAAAAAAAAAAAA  A  A  A  A  A  4 

4  4  4  AAAAAaaAAAAAAAAA  4  4  4  4  4 

4  4  4  4AAAAAAAAAAAA4A  A  4  4  4  4  4 

A  A  A  AAAAAAAAAAAAAA4  4  4  A  A  4  4 

A  A  A  AAAAAAAAAAAAA44  A  4  4  A  4  4 

A  AAA  ^AAAAAAAAAAAff  44  4  4  4  4 
A  A  A  AA A A A AAAA44444 4  4  4  A  4  4  4 
A  A  A  A AAAAAAAAAAAA4A4  A  A  A  4  4 
A  A  A  A  AAAAAA4444444  4  4  4  4  A  4  4 

A  A  A  A  A  AAA  AA4A444444  4  4  4  4  4  4 


A  A  AAAAAAA444444A4 4 4  4  4  A  4| 


A  A  A AAAAAAA4444444 4  4  4  4  4  4 


§  s 
fc  I 

H  " 

§  f 


A  A  A4AAA4A44444444  4  4  4  4  4  4] 


AAAAAaAAAAA  A  A  a  a 

AAAAAAAAAAAA  a  a  a 

aaaaa  aaaaaa  a  a  a  a 

AAAAAAAAAAA  A  A  A  A 

aaaaaaaaaaaat  r  r 
aaaaaaaaaaa  r  r  t  t 

aaaaaaaaaa  t  t  T  f  t 

aaaaaaaa  r  A  T  t  <  <  < 

aaaaaaaa  a  a  t  <  <  «•  1 


*  A  A  A  A  A  A  A  A 
A  A  A  A  A  A  A  A  A 
A  A  A  A  A  A  A  A  A 
A  A  A  A  A  A  A  A  A 
A  A  A  A  A  A  A  A  A 
A  A  A  A  A  A  A  A  A 


77  zJ j  4 

Ut .  *  > 

It.**  , 


J; 

,  ;TW 

: 

: : . . .  '>v^aaaaa 

<■■  -  sS,SS,Ss 

,  s*  A,  A  A  A 

JJT~- 

^  ^  <9  A  4  4  4  * 


A  4  4  4  «4  y 

4444444  4 


4  4  4  4 


4  4  4 


44444444 

44444444 


4  4  4  4 


4  4  4 


4444444  4 

4444444  4 

4444444  4 


A  A  A44444A4444**** '4  Y  4  , 
A  A  4  444444444+44««  *  *  4  t  i 
A  ♦  f  f  *  i 

A  A  ♦  ♦  ttf  ♦***#*»#*«  «  ,  ,  ,  , 

A  ^  tttf  f  f  f  A  T  ' 

T  *  If  , 

9  1  kkt^^/^*******  «  V  V 

f  I  J  t£t  •  9  *  9 

"T  / /XZ tvv^mmmi  i  .»  a 
/  . . 

»  •  4 

>1  «  «  4 

^  ^9999^999  44  4 

▼  t  ^99999994  4  4  4  4 

<  ««t>999V49V994  4  4  4  4 
4  4«.l<.9*kvvmw«4  4  4  4  4 
i.  A  lLU4>l<vvyvYWVV  4  4  4  4 
A  A  A  A  A  AAWw  WVW4  4  4  4  4  4 
A  A  AAAAAVVVVW4VV4  4  4  4  4 
y  V  VVVVWVVWV W44  4  4  4  4 
V  y  VWVVVVVVVV444  4  4  4  4  4 
y  v  yyyyyy yyyyyyyyy  4  4  4  4 
4  4  44444yyyyy444444  4  4  4 

4  4  4444444444444444  4  4  4  ' 

4  4  4444444444444444 4 4  4 

4  4  4444444444444444  4  4  4 

4  4  4444444444444444  4  4  4 

4  4  4444444444444444  4  4  4 

4  4  4444444444444444  4  4  4 


4  4  444444444444444  4  4  4 


4  4  444444444444444 4 4  4 


4  4  4444444444444444  4  4  4  4 


4  4  4444444444444444  4  4  4  4 


444444  4  444  44444444444444  4  4  4  4  4  4 


4  4  4  4 


4  4  4  4  4  4444444444444444  4  4  4  4  4 


^W«'<j55SiI5i 


MTV)  Mft  minOVA 


vvvvvvvvv 
vwvyvvvvv 

t 

*  V  V  V  V  V  V  V  V 

v  *  v  v  v  v  v  v  y 

v  y  v  v  v  y  y  y  v 

y  y  w  v  y  v  y  y  y 

yyvyyyy  y  y 
yyyyyyyyyy^ 


v  V  yy^^k-wUUUUt  4 

*  Vk»y»^^UUi4lil  4 

V  4  i 

4  y  «v»rv»»»»t4AA14i4  4  4 

4  «  ♦»»♦»»»>*>«  »44 4  4  4  4  4 

1  *  «***«»•““““  4  4 

*  yw»««»****“**“ 


vvyyyy  y  y  v 
vvyyyyyyy 
vvyyyyyyy 
y  w  l.  i.  4.  y  y  y  y 

y  y  i.  4,  y  y  y  y  y 

y  4.  y  A-  y  y  y  y 


y  y  y  y  y 
y  y  y  v  t 
y  v  v  y  y 


y  4j^ 

*  [7"> 

r5-^  / 

-V'  ✓  ^ 

//  ✓ 
S  /  ' 


r= ^ 
y 

✓  / 


r-1' 

/  / 
/  / 
/  / 
✓  ■" 

',A 


JLlUUW**-**** 
J  •  ///W/********* 
*  /////^‘it4>  ***"'** 
s  SAfAM********* 


L?T«if<ij.um>>  i  »  » 
[7 ^IWUMMt  »  t 
. . . 


4  »  » 


r /'S/S/S  ' 

y  y^y  /  ^  V 

VS//  /  '  A*  1 

VY/  /  /  4 

yy  /  /  4  4 

/yv  /  /  »  fc 

/yy  /  ju  *•  r  *  -« 

rj  /  *  -  ‘  T  r  “  a  * 

///•*'*'■*** 

4  1/LlV  A  *  *  - 


*  A  A  A  A 


r*“***‘**“  *■ 
>  >  »»11**«“**‘»»»  > 

»  »  ►  * 

7  >  »»>»T 777<««**A*»»  »  ► 

t  •»  17T7l1^144AAAAA*  4  4  4 
«f  1  1117l1AyiAMAA4  4  4  4  4 
A  A  i7A111AAAAAAA*4  4  4  4  4 
A  A  AAAAAAAAAAAAMA4  4  4  4 
A  A  AAAAAAAAAAAAA4A4  4  4  4 
A  A  AAAAAAAAA AAAAA44 444 
A  4  AAAAMAAAAAAAAA4  4  4  4 
A  A  AAAAAAAAAAAAA4A4 444 
A  A  A  AAAAAAAAAAAA4  4  4  4  4 
A  A  AAAAAAAAAAAAAA4  4  *  *  * 


A  A  AAAAAAAAA4A4A44  4  4  4  4 


A  A  AAA4A4AA4444444  4  4  4  4  4  4 


A  AAAAAAA44444444  4  4  4  4  4  4 


A  A  A AAAAAAAAAAA^^  41  4  4  4  4 


A  A  A44A44444444444  4  4  4  4  4  4 


A  A  A  444444444444^444  4  4  4  4  4 


VMUINA  MY  RUN 


44444a  A  A  A  A  A  A  A 

aaaaaaaa  a  a  a  a  a 

A  A  A  A  A  AAA  A  A  A  A  A 


aaaaaaaa  a  a  a  a  a 

AAAAAAAA  A  A  A  A  A 

aaaaaaaa  a  a  a  r  r 

44444AAAArrr  r 
4  4  4  4  a  a  a  a  A  r  t  *  < 

4  4  4  4  4  A  A  a  a  T  X  <  * 

4  A  A  A  A  A  A  A  4  T  «  < 


1  1  ^  1  ^ 

y  a  a  a 


y  a  a  a 
y  y  a  a 


’  i>  0 


/////f\  /Tr  ■•>»>»»  j 

77V  >>>>»* 

/  / '  '  *  ■*  >>>>*■» 

'f  '  '  *  «  >>>»*•» 

'  ^  <  1  «  f  >  >  >  >  t  J 


A  A  4AA44441M444414  4  A  A  A  A 
A  *  44444444444444111  A  A  A  A 
f  4**44444444«ii1f41i,, 

4  ♦  t  ttff#4444Ml4li  i  4  i  i  i 
r  %  , 
,  ^  tttff 

T  ^  I  f  ttffff0*m9wm  W  •  ¥  V  * 

■  f  m  a  »  »  i 

/  «  »  4 

1114  *  4 

/»/>!wm4^1||  n  t  ,  ,  ,  | 


/Via 


>>>444 


'H  *  *  «  4  m  m  •  4  4  4  A 

I  *Wfc**t”**"  •  4  4  * 
•  4*+«M*4*«  4  4  4  4  4 
*♦♦*411414  4  4  4  4  4  4 

. . 4 

4  <  <<<(W*v?VVVTtin  4  4  4  4  4 

4  4  4  4  4  4 

4. 

*•  *•  ttttH-yywyvyyV  4  4  4  4  4  4 

A-  >  >AUvwvwyy»v4  4  4  4  4  4  4 

V  V  WVW>VWVWVW44  4  4  4  4  4  4 

V  V  VWyyvVVWVV444  4  4  4  4  4  4  4 

V  V  vvwm VWVW4V4  4  4  4  4  4 

y  V  VVVVVVVVVVVV44  4  4  4  4  4  4  4 
y  V  V4444V44444V444 4 4  4  4  4  4 
4  V  444444444444444  4  4  4  4  4  4 
4  4  4444444444444444  4  4  4  4  4 

4  4  444444444444444  4  4  4  4  4  4 

4  4  4444444444444444  4  4  4  4  4 

4  4  4444444444444444  4  4  4  4  4 

4  4  444444444444444  4  4  4  4  4  4 

4  4  4444444444444444  4  4  4  4  4 


4  A  A  4  4  4  4 


4  4  4444444444444444  4  4  4  4  4 


4  14  4  4 


4  4  4  44444444444444  4  4  4  4  4  4 
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Plate  22.  Difference  between  existing  conditions  and  alternative  d;  peak  flood 
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Plate  28.  Wave-induced  current  pattern 
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Plate  37.  Superposition  of  tidal  and  wave-induced  currents  at  peak  flood 

H  =  10  ft,  T  =  II  sec,  6  =  WSW 
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Plate  50A.  Superposition  of  tidal  and  wave-induced  current  patterns  ai  peak  ebb 
for  alternative  b:  H  =14  ft,  T  =  11  sec,  0  =  WNW 
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